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1 FORWARD ERROR CORRECTOR 

CROSS-REFERENCE TO RELATED APPLICATION 

This patent application claims the benefit of the filing date of U.S. Provisional Patent 
5 Application Serial No. 60/107,879, filed November 9, 1998 and entitled FORWARD ERROR 
CORRECTOR, the entire contents of which are hereby expressly incorporated by reference. 

BACKGROUND OF THE INVENTION 
1 . Field of the Invention 
10 The present invention relates to an apparatus for correcting errors present in stored or 

transmitted data; and, more particularly, to an apparatus for evaluating an error evaluator 
polynomial, an error locator polynomial and a differential polynomial which are used in 
correcting errors in the data encoded by using an algebraic code, such as a Reed-Solomon code. 

1 5 2. Description of Related Art 

Noise occurring during a process of transmitting, storing or retrieving data can in turn 
cause errors in the transmitted, stored or retrieved data. Accordingly, various encoding 
techniques, having the capability of rectifying such errors, for encoding the data to be transmitted 
or stored have been developed. 

20 In such encoding techniques, a set of check bits is appended to a group of message or 

information bits to form a codeword. The check bits, which are determined by an encoder, are 
used to detect and correct the errors. In this regard, the encoder essentially treats the bits 
comprising the message bits as coefficients of a binary message polynomial and derives the check 
bits by multiplying the message polynomial R(x) with a code generator polynomial G(x) or 

25 dividing R(x) by G(x), to thereby provide a codeword polynomial C(x). The code generator 
polynomial is selected to impart desired properties to a codeword upon which it operates so that 
the codeword will belong to a particular class of error-correcting binary group codes (see, e.g., 
S. Lin et al., "Error Control Coding: Fundamentals and Applications", Prentice-Hall, 1983). 
One class of error correcting codes is the well-known BCH 

30 (Bose-Chaudhuri-Hocquenghen) codes, which include the Reed-Solomon ("RS") code. The 
mathematical basis of the RS code is explained in, e.g., the aforementioned reference by Lin et 
al. and also in Berlekamp, "Algebraic Coding Theory", McGraw-Hill, 1968, which is further 
referred to in U.S. Pat. No. 4, 162,480 issued to Berlekamp. The aforementioned references are 
hereby incorporated by reference in pertinent part. 

35 

SUMMARY OF THE INVENTION 

The invention herein provides a method and apparatus for decoding an algebraic-coded 
message. The method can include the steps of determining a discrepancy indicator, with the 
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discrepancy being between a calculated and a predicted value; determining an error locator 
polynomial using a selected class of error correction algorithms, such as, for example, a 
Berlekamp-Massey algorithm; and detecting an uncorrectable message using the selected error 
correction algorithm. The apparatus is composed of storage devices which can include 
recirculating storage devices; arithmetic components attached to the storage devices, the 
components operating over a Galois Field on selected contents of the storage devices; and an 
uncorrectable message detector, connected with the storage devices and the arithmetic 
components. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is an illustration of a algebraic decoder according to the invention herein; 

Figure 2 is a data flow diagram of a modified Berlekamp-Massey algorithm according to 
the invention herein; 

Figure 3 is a block diagram illustrative of an exemplary embodiment of the present 
invention; 

Figure 4 is a block diagram of a circular syndrome generator according to the present 
invention; and 

Figure 5 is a block logic diagram of a logic register which can be used in the circular 
syndrome generator illustrated in Figure 4. 

EXEMPLARY EMBODIMENTS OF THE INVENTION 

The invention herein provides an apparatus for and a method of decoding algebraic codes, 
including BCH codes, and more specifically, Reed-Solomon codes, such that uncorrectable 
messages, or portions of received encoded data, are detected. Furthermore, the invention herein 
provides for a more area-efficient device implementation of the aforementioned method. For the 
purposes of illustration, the present invention will be described in terms of a subset of the BCH 
codes, namely Reed-Solomon (RS) codes. 

The Reed Solomon (RS) encoding technique appends to each block of k user data 
symbols 2t redundancy symbols to create an encoded message block (where / represents the 
designed symbol error correcting capacity of the code). These 2t symbols, or elements, are 
selected from the Galois Field to be the roots of the generator polynomial. Therefore, there are 
jfc+2/ symbols in a RS-encoded message block. The entire message block is viewed as a 
polynomial and evaluated as a polynomial at some Galois Field element. The Galois Field 
element at which the polynomial is evaluated will be located at one roots of the generator 
polynomial that are used to create the RS code. The RS code views the w-bit symbols as elements 
of a Galois Field (GF(2 n )). A Galois field is a finite field, the elements of which may be 
represented as polynomials in a, where a is a root of an irreducible polynomial of degree n. The 
RS codeword consists of a block of /7-bit symbols. Typically, n = 8 and the 8-bit symbols are 
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1 referred to as bytes. Constructing the Galois field GF(2") requires a defining polynomial F(x) of 
degree n. In addition, a primitive element ft is chosen so that every nonzero element of GF(2") 
is a power of fi. The element /?is not necessarily a root of F(x). 

A RS codeword C is viewed as a polynomial C(x) and the redundancy symbols are chosen 
5 so that the roots of C(x) include the roots of a generator polynomial G(x) whose roots are 2t 
consecutive powers of p. The k user data symbols are viewed as the high order coefficients of 
a degree k+2t-J polynomial, and the redundancy symbols are the coefficients of the remainder 
when this polynomial is divided by G(x). 

The process of corrupting the original code block C(x) with errors can be viewed as 
10 adding an error polynomial E(x) to C(x). The resultant corrupted polynomial is known as the 
received polynomial R (x), where R(x)=C(x)+E(x). The v non-zero terms of the error polynomial 
contain the necessary information to completely reconstruct the original data C(x), since each 
term corresponds to a symbol error location and magnitude. 

Typically, RS decoding is a tripartite analysis: (1) syndrome computation; (2) solution 
1 5 of the error magnitude and locator polynomials; and (3) error location and magnitude estimation 
by respective implementations of, for example, a Chien search and the Forney algorithm. The 
syndromes contain error information divorced form the actual information that is intended to be 
analyzed for errors. The error locator polynomial provides information regarding the location 
of an error in the received signal, and the magnitude of the error can be determined by using both 
20 the magnitude and the locator polynomials. 

The thrust of the RS error correction procedure is to reconstruct the error polynomial E(x). 
Three polynomials are used to correct a received polynomial R(x): S(x), a syndrome polynomial; 
A(x), an error locator (or error location) polynomial; and M(x) an error magnitude polynomial. 
The syndromes are computed by evaluating the polynomial R(x) at all roots of G(x). These values 
25 are called syndromes and the syndrome polynomial S(x) has these values as coefficients. The 
syndrome polynomial S(x) is used to determine the existence of errors. The error locator 
polynomial A(x) and the error magnitude polynomial M(x) are computed from S(x) by a key 
equation solver. The roots of the error locator polynomial A(x) indicate positions in the data that 
are erroneous and both the error locator polynomial A(x) and the error magnitude polynomial 
30 M(x) are used to determine the true values of the erroneous data. 

Two frequently-used RS error correction algorithms are the Berlekamp-Massey and the 
Euclid algorithms. The present invention recasts the Berlekamp-Massey algorithm such that the 
inversion process typically associated with the traditional Berlekamp-Massey (tBM) algorithm 
is eliminated. This is important because the inversion process includes determining the 
35 reciprocal of certain Galois field elements using division. Division is a time consuming 
arithmetic operation, the implementation of which can occupy needed component area in a device 
design. Therefore, the present invention can be particularly advantageous where area-efficient 
layout of a decoder device is desirable. 
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1 For further elaboration of the decoding process over Galois fields, including tBM, Chien 

searching, and Forney's Algorithm, see Theory and Practice of Error Control Codes by Richard 
E. Blahut (Addison- Wesley, 1983) which is incorporated by reference in pertinent part herein. 
Figure 1 illustrates an implementation of this algorithm, in which a raw received signal 

5 1 is directed to RS decoder unit 2 that is used to determine the error locations and error values. 
Signal 1 is provided to syndrome generator 3 and delay unit 4. In syndrome generator 3, the 
several syndromes 5 associated with the selected encoding are derived and transmitted to 
polynomial solver 5. The syndrome generator 3 calculates one syndrome for each of the 2t roots 
of G(x). Polynomial solver 6 utilizes the syndromes to determine the coefficients of the error 

1 0 location polynomial A(x) 7 and the coefficients of the error magnitude polynomial M(x) 8, which 
in turn are transmitted to error estimator 9. Estimator 9 calculates error signal 10 which is 
combined in summer 11 with delayed raw received input 12 to provide corrected data 13. 
Estimator 9 can include Chien search unit 14 which utilizes the error location polynomial A(x) 
to search for the roots of the error locator polynomial, r h . . . , r tf . Typically, the Chien search unit 

15 14 employs a root finding technique which involves evaluating the error locator polynomial at 
all elements in the field GF(2 n ). The roots of the error locator polynomial r„ . . . , r v determine 
the error locations. The error values are then determined using Forney's algorithm unit 15. The 
delayed raw received input 12 is then corrected using the output of the Forney algorithm unit 15 
and the raw received input which is transmitted by delay unit 4. 

20 Traditionally, the Berlekamp-Massey (tBM) algorithm, which usually is realized in 

polynomial solver 6 can described by: 

25 

lW,(r-4-, + (l-* r )4-i (2) 



~A (0 (x)~ 




' 1 -A,x 


>-'>(x) 


_B {r) (x) 




_A"X (1 -<*>)*_ 


_B (r - l) (x) 



r=l,..., 2t where 6 r = 1 if both A, * 0 and 2L,_, <. r-1, and otherwise 6 = 0. Then A m (x)\s the 
smallest-degree polynomial with the properties that Af' )= \, and 

S r + X A > ' )S r-J = 0 r = L 2t + 1 >-> 2t 
7=1 

where initial conditions are A (0> (x)-\, BP°fr)=l, and L o =0. 

It is evident that the inversion indicated in Eq. 3 requires a division operation. 
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1 The tBM algorithm is capable of properly decoding messages that can be decoded 

properly, however if there is an uncorrectable case which is detectable as being uncorrectable, 
the uncorrectable error may be undetected and the message decoded as if it did contain a 
correctable error. Many times, this improperly decoded message can create additional difficulties 

5 because the error may propagate through other processes in the system which employs tBM. 

According to the invention herein, the modified Berlekamp-Massey (mBM) can be 
described by the following equations: 

10 

A,. = A_A,._, + x A (5) 



15 



20 



B, = \ 1 ( 6a ) 

(6b) 



where: A * 0 

5 0 =1 
4^1 



Utilization of mBM for RS decoding can be advantageous because: (1) inversion is 
eliminated; (2) the control structure associated with the mBM algorithm is simplified relative 
to that of tBM; and (3) the termination conditions of tBM are modified such that if the code is 
uncorrectable, errors otherwise undetected by tBM, are detected and flagged as such. 
25 One implementation of the mBM algorithm is as follows, as represented in Pascal code: 



PROCEDURE FindLocatorBMC( VAR 

VAR Cnt: O..MaxParitySyms-l; 
Pwr: CMaxParitySyms; 
State: (Alpha,Beta); 
Deg: INTEGER; 

Del, DelO: Words; 
J: INTEGER; 

TempPoly: Polynomial; 
B : Polynomial; 

BEGIN 
B.L. :=0;B.D [0] :=1; 
Locator.L := 0; Locator.D [0] :=1; 
Deg := 0;Pwr :=0; DelO := 1; 
State := Alpha; 



:,Locator:Polynomial; VAR OK.BOOLEAN ); 

{ Loop Index 
{ Power Counter 
{ State Machine State 
{ Degree 

{ Discrepancies 
{ Del Index 

{ Temporary Polynomial 
{ Secondary Polynomial 

{ BEGIN FindLocatorBMC 
{ Initial B 

{ Initial Locator Poly 
{ Cntr Initialization 
{ Machine State 
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FOR Cnt := ParitySyms-1 DOWNTO 0 DO BEGIN 



10 



15 



20 



{ Algorithm Loop 
{ Calculate Del 



Del := 0; 

FOR J := 0 TO LOCATORL DO 
IF SyndromeX >= (ParitySyms-l-Cnt-J) THEN 

Del:= Add( Del, Multiply( Locator.D[J],Syndrome.D[ParitySyms-l-Cnt-J])); 

TempPoly := { Do Common Update 

PolyAdd( WordTimes( Locator, DelO ), PolyShift( WordTimes( B, Del ), 1) ); 



25 



IF (State=Alpha) AND (DeloO) THEN BEGIN 
{ writeln( stderr, ' B<-L');} 
B := Locator; 
DelO := Del 
END 

ELSE BEGIN 
{ writeln( stderr, ' B<-xB' );} 
B:= PolyShift(B, 1) 
END; 

IF State=Alpha THEN BEGIN 
IF Del=0 THEN Pwr := Pwr +1 
ELSE State := Beta 
END 

ELSE BEGIN 
Deg := Deg+1; 

IF Pwr =0 THEN State := Alpha 

ELSE Pwr := Pwr-1 

END; 

Locator := TempPoly 
END; 

Locator := PolyDenormalize( Locator, Deg); 

OK := State=Alpha 

END; 



{ Do Step A 

{ Do Step A 
{ Do Step B 

{ Do Step B 

{ State is Alpha 
{ Increment Power Cntr} 
{ Update Next State 
{ State is Alpha 

{ State is Beta 

{ Update Next State 
{ Decrement Power Cntr 
{ State is Beta 

{ Update Locator 
{ Algorithm Loop 

{Update Locator Degree 

{ END FindLocatorBMC 



Often, when a forward error corrector properly detects an uncorrectable error, the 
existence of such an error usually is verified in a process by which the syndrome polynomials 

30 are recomputed. This approach can carry a substantial penalty relative to the process efficiency. 
Instead, an embodiment of the invention herein, having an improved control structure, verifies 
the existence of an uncorrectable error by checking the state of polynomial solver 6 at the end 
of the polynomial solving process. 

Figure 2 exemplifies an embodiment of the process 20 implementing the aforementioned 

35 improved control structure in the context of the mBM algorithm recited, in Equations 4, 5, and 
6(a)-(b). Although the implementations described herein are postured for standard RS codes 
having a block length of, for example, 255 elements, such implementations also may be used 
in the context of extended Reed-Solomon codes which, in the example herein, would have 256 
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1 elements in the message block, i.e., have 256 elements in associated the Galois Field. It is 
desirable that, in step 2 1 , the control variables DEG, P WR, and STATE, as well as error locator 
variables be initialized. It further is desirable to iterate through steps 23, 24, 25, and 26, 2t 
times, where 2t is the number of syndrome polynomials to be evaluated, and t is the error 

5 correcting capability of the preselected code. Thus, at step 22, a counter tracking the number 
of completed iterations is employed. No additions or subtractions are needed in implementing 
the control variables, and only count up or down functions are used. Step 23 essentially 
implements Equation 4, in which the discrepancy value DEL, associated with a particular 
iteration, is determined. Similarly, step 24 implements Equation 5 in which the error locator 

1 0 polynomial is updated. In step 25, auxiliary polynomial is updated according to Equation 6a 
in substep 27, or Equation 6b in substep 28, based on conditions determined by logic 26. For 
logic 29, it is desirable for both STATE = ALPHA AND DEL < > zero to direct the data flow 
via an implementation of Equation 6a in substep 27; otherwise substep 28 is used, 
implementing Equation 6b. Unlike the tBM algorithm where the polynomial shift term 

15 ) x in Equation 3 has been normalized, the mBM algorithm does not require 

normalization, avoiding an inversion/division operation. After the auxiliary polynomial is 
updated in step 25, the controller state is updated in step 26. 

In general, the degree of the error locator polynomial is tracked by DEG, which is an 
upcounter descriptive of the true degree of the error locator polynomial and, thus, the number 

20 of errors in the message block. It also is desirable to construct an error locator polynomial who 
roots equate to the locations of an error. Essentially, process 20 attempts to synthesize a linear 
feedback shift register (LFSR) that predicts the values of the syndrome polynomial. Such a 
LFSR can be useful to find the error locations. Discrepancy value, DEL, then exposes a 
discrepancy between the predicted value of the syndrome polynomial, and the value of the 

25 current syndrome polynomial, and invites further processing to discover the location of the 
errors. PWR is a counter that keeps track of the number of times that the controller previously 
remained in STATE = ALPHA. It is desirable to have the STATE remain in control state 
BETA for a count equivalent to the number of times that STATE previously remained in control 
state ALPHA. 

30 For the purposes of the invention herein, STATE can be used to (1 ) determine whether 

the error correction analysis ought to follow the flow of Equation 6a or 6b; (2) assist in 
determining whether the degree of the error locator polynomial ought to be increased; and (3) 
whether there exists an uncorrectable error. At the end of 2t iterations, the value of STATE is 
once again determined at step 35 . If the result is STATE = ALPHA, then the code is potentially 

35 valid; on the other hand, if STATE = BETA, then the error is flagged as uncorrectable. 
Potentially valid codes where STATE = ALPHA at step 35, also can be subjected to additional 
validation before being subsequently decoded. Indeed, in one subsequent operation, the number 
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1 of the error locator polynomial zeroes is compared with the value of DEG. A discrepancy 
between these two values also is indicative of an uncorrectable error. 

Figure 3 is an exemplary embodiment of a polynomial solver using the mBM algorithm. 
Solver 50 can include syndrome generator register 5 1 , recirculating syndrome register 52, first 
5 delay register 53, locator polynomial^, Register 54, auxiliary polynomial (5,) register 55, 
seconddelayregister56,firstmultiplier57,secondmultiplier58,adder59, A register 60, and A_ 
register 61. In another embodiment, syndrome generator 51 can be separate from solver 50. 
It is desirable for multipliers 57, 58, and adder 59 to operate on Galois Field elements. It also 
is desirable for register 51 to be logically arranged as a circular register or loop, such that 
1 0 particular register values can be used in a pre-defined sequence. Furthermore, it is desirable that 
registers 52, 54, and 55 be logically arranged as push-down stacks or FIFOs, and also that the 
values contained therein rotate synchronously. The error locator polynomial A, are arranged 
in register 54 such that the least significant coefficient is at the top and the stack "grows down" 
as subsequent values are determined. It is desirable for the syndrome recirculating register 52 
1 5 to operate such that the least significant coefficient is aligned with the bottom of the register, 
and the most significant with the top. 

In the example of Figure 3, the error correcting capability, t, is selected to be 5 elements 
and, thus, syndrome register 51 is designed to employ 2t, or 10, individual elements. 
Additionally, register 52 is selected to use / elements, register 54 is chosen to employ /+/ 
20 elements, and register 55 is intended to operate with / elements. 

Initially, recirculating syndrome register 52 is pre-loaded with zeroes. As expressed in 
the aforementioned algorithm, an initial step involves calculating the discrepancy value DEL, 
which can be stored in register 60. A previous value for DEL, DEL0, is provided in register 6 1 . 
To calculate the initial value for DEL, first syndrome value S 0 is shifted into register 52, which 
25 value is received in first multiplier 57 along with the initial value of DEL, namely, DEL0,and 
the Mh value in locator polynomial register 54. After the indicated multiplication and creation 
of a DEL value, the values in registers 52, 54, and 55 are shifted down by one element. At first, 
the values in register 52 are zero, however, with subsequent clocking, successive values of 5, 
enter register 52 and are recirculated therethrough, for the clock cycles equivalent to /' = 0 to 21- 
30 1. As S 0 exits register 52 into first multiplier 57, corresponding values of A 0 are also 
transmitted to first multiplier 57 such that the value S 0 A 0 is determined. 

Concurrently with this calculation, value B 0 from register 55 is multiplied with then 
extant value for DEL in second multiplier 58 and the result is summed with S 0 A 0 in adder 59 
to produce the next value for DEL. This value of DEL is used to produce the next value for the 
35 error locator polynomial, namely, A,. After this calculation, value S 0 is recirculated such that 
it bypasses first delay register 53 and re-enters recirculating register 52 at the top of the stack 
during the next clock cycle. During this next clock cycle, syndrome value S, is aligned, and 
multiplied, with A 0 , and S 0 is aligned, and multiplied, with A, . This process repeats such that 
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1 each of the syndrome values properly iterates through in the evaluation of the locator 
polynomial. 

With the above information, a skilled artisan would be able to see the manner in which 
the values for error locator polynomial A, and auxiliary polynomial B, are determined. Where 
5 it is desired to rotate the values of /?, through register 55, second delay register 56 is bypassed. 
On the other hand, where it is desirable to utilize a previously calculated value of B i9 as 
indicated in the aforementioned algorithm, the value of £, is directed into second delay register 
56. 

Continuing in Figure 3, calculation of the magnitude polynomial will be described. At 

10 the completion of 2/ iterations as described above, register 52 will contain values S 4 -S 8 In 
essence, determination of the magnitude polynomial can be modeled as 
M(x) ~ A(x) S(x) mod JC 2 ', in which the multiplication will be truncated after the 2t term. 
Indeed, only t terms need be determined under the assumption that no uncorrectable error was 
encountered. During the final iterations of the calculation of the error locator polynomial, 

1 5 register 52 is loaded with zeros such that, at the completion of 2t iterations, all storage locations 
in register 52 contain a zero value. After 2t iterations, the values in register 5 1 will be restored 
to their original positions, register 52 will contain all zero values and register 54 will contain 
the error locator polynomial, A, . In am manner similar to the computation of the locator 
polynomial coefficients, the error magnitude coefficients are calculated iteratively. After / 

20 iterations, S 0 will be at the logical bottom of register 52 and A 0 at the logical bottom of register 
54. At the completion of /+/ iterations, the product of multiplier 57 will be S 0 A 0 , the first 
term of the error magnitude polynomial. The output of adder 59 is directed to the logical top 
of register 55, which now will be used to build the magnitude polynomial. After iteration t+2, 
syndrome value S 0 will be aligned with locator value A, , giving the product S 0 A, ; syndrome 

25 value S 1 will be aligned with A 0 , giving the product S, A 0 ; the summation of which giving 
S 0 A 1 + S f A 0 , which is the second term of the error magnitude polynomial. This process will 
continue until all values of M(x) are so determined. At iteration 2/, all of the coefficients for 
the error magnitude polynomial will have been calculated. At this point, data flow of the error 
locator polynomial in register 54 and the error magnitude polynomial in register 55 can be 

30 directed out of solver 50. 

Figure 4 illustrates one embodiment of a circular syndrome generator 70 that can be 
employed as register 5 1 in Figure 3, modified to accommodate an error correcting capability of 
/=$. Figure 5 is an embodiment 72 of one of the individual registers 71 in circular syndrome 
generator 70 in Figure 4. Although a circular syndrome generator is shown, it is by no means 

35 the only form of syndrome generator that can be employed as register 5 1 . 

The foregoing merely illustrates the principles, of the invention, and it will thus be 
appreciated that those skilled in the art will be able to devise various alternative arrangements 
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What is claimed is: 

1 . A method for decoding an algebraic-coded message, comprising the steps of: 

a. determining a discrepancy indicator; 

b. determining an error locator polynomial according to an algorithm from a 
selected class of decoding algorithms wherein an uncorrectable message is detected; and 

c. producing a perceptible indication of the detected uncorrectable message. 

2. The method of claim 1 wherein the selected class of decoding algorithms is a 
Berlekamp-Massey algorithm. 

3 . The method of claim 2 wherein the Berlekamp-Massey algorithm is an iterative 
algorithm. 

4. The method of claim 3 wherein the iterative algorithm further includes the step 
of evaluating the error locator polynomial for the presence of an uncorrectable message after 
an iteration. 

5. An apparatus for decoding an algebraic-coded message, comprising: 

a. a plurality of storage devices, selected ones of the storage devices being 
recirculating storage devices; 

b. a plurality of arithmetic components operably connected with the plurality of 
storage locations, the arithmetic components being operative within a Galois field; and 

c. an uncorrectable message detector, operably connected with the storage devices 
and the arithmetic components. 

6. The apparatus of claim 5 wherein one of the recirculating storage devices is one 
of a recirculating syndrome polynomial storage device, a recirculating error location polynomial 
storage device, and an auxiliary coefficient storage device. 

7. The apparatus of claim 5, further comprising a syndrome generator. 

8. The apparatus of claim 7 wherein the syndrome generator is a circular syndrome 
generator. 

9. The apparatus of claim 8 wherein the circular syndrome generator is a 
recirculating circular syndrome generator. 
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1 0. The apparatus of claim 6 wherein the one of a recirculating syndrome polynomial 
storage device, a recirculating error location polynomial storage device, and an auxiliary 
coefficient storage device, is an area-efficient device. 

11. An apparatus for decoding an algebraic-coded message, comprising a circular 
syndrome generator. 
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